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In this paper we present a new simple method of construction of infinite number of solutions 
of Freidmann equations from the already known ones, which allows for a startling conclusion of 
practical impossibility of correct predictions on the universe's future dynamics which are based 
solely on astronomical observations on the value of a scale factor. In addition, we present particular 
examples of newly constructed solutions, such as the ones, describing the smooth dynamical (de)- 
phantomization, and the models lacking the events horizons (both in classical and brane world 
cases) . The generalization of the method to the simplest anisotropic universes are presented as well. 
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I. INTRODUCTION 
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Imagine for a second an immortal astronomer, rest- 
lessly surveying an expansion of the universe, how it goes 
during 10 20 consequent years. Suppose that she already 
knows for sure about the flatness of her universe. Then, 
as time goes, the astronomer will gather an (impressive) 
amount of observational data, visually presented as a set 
of points on the graph {U, a(U)}. If those points are lo- 
cated with appropriate density, then it should be possible 
to choose some continuous function a(t), which with high 
precision will approximate the gained set. Of course, this 
function must not be unique, but lets assume for now 
that our astronomer have managed (not without some 
efforts) to find out the most satisfactory one; when this 
is done she might conclude that now she posses the ac- 
tual time dependence of observed scale factor. Moreover, 
since both matters density and pressure are expressed in 
terms of scale factor (and its derivatives) explicitly (we 
remind here, that under the assumption the sign of cur- 
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then she might as well conclude that the observations 
have given her the values of all basic characteristics of 
the universe. Can we say now that the further evolution 
of her universe will be completely defined? No, we can't. 

This paper contains the discussion of an extremely 
simple method which allows one to construct new so- 
lutions of the cosmological Friedmann equations starting 
out from the previously known (and more simple) ones. 
The method itself will be denoted as the method of lin- 
earization, and the transformations which allows one to 
construct new solutions starting out from another ones 
will be refereed to as "the dressing" or as a "dressing 
procedure". What is new and interesting about such 
transformations is that they admit the invariants. More 
precisely, the special linear combination of the density p 
and pressure p (of a form U = ap + (3p) is covariant with 
respect to these transformations {ap + ftp — ap + /3p). 
(For example, in one particular case p = p). 

There are in total four aims of this paper. First of 
all, we will show that for any function a(t) there exist 
such solutions a n , which, being similar to a(t) at given 
time interval with any given precision rate, will com- 
pletely differ from a{t) later (note here, that possibility 
of such scenario, i.e. of indeterminacy of universe fu- 
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ture has been previously noted by A. A. Starobinsky in 
These solutions will be constructed via the dress- 
ing procedure starting out from the initial solution a(t) 
(Sec. II). Secondly, as we shall see in Section III, upon 
further usage of a dressing procedure one can construct 
special solutions, describing the smooth dynamical (de)- 
phantomization. Thirdly, in Sec. IV we will show that 
described transformations allow one to construct space- 
times without event horizons. Finally, in Sec. V we will 
take a brief look on possible generalizations of presented 
method for the simplest anisotropic cosmological model 
and the particular Randal-Sundrum-I brane model and 
will show that in both cases the similar dressing proce- 
dures would still be valid and they yield the same results 
as the ones described above. 



mathematical methods which were developed for the lin- 
ear differential equations. This is the reason why we call 
our approach the method of linearization. The crucial 
point of this paper is connected to the simple general- 
ization of results above. More precisely, the following 
proposition is hold: 

Proposition. Let a = a(t) (with p — p(t), p = p(t)) 
be the solution of ^ with k — 0. Then the function 
ijj n = a n is the solution of the Schrodinger equation 



with potential 



i — u rn 



3n 



U n = n p - — (p + p) 



(6) 
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II. DRESSING PROCEDURE 



Let us write the Friedmann equations as 

2 
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where k = 0, ±1. Throughout the paper (excluding 
Sec.V) we'll stick to the metric units with 8nG/3 = c = 1. 
If the universe is filled with a self-acting and minimally 
coupled scalar field with Lagrangian 



then the energy density and pressure are 



(3) 



p = K + V, p = K-V, 



therefore 



V=\{ P -P), K=\(p+p). (4) 

Our starting point is that the volume function ip = a 3 
satisfies a simple second-order differential equation ((2j) 



(5) 



In ([5]) the potential V is represented as a function of 
time t. For simple forms of the potential one can find 
the general solution of ||5j), containing both the solution 
used for the construction of this potential and a lineary 
independent one as well. Substituting this general solu- 
tion into the @ one can calculate p and p. Then using 
(f4|) one gets the new potential V such that V(t) = V(t) 
but whose form is different from V if V and V are rep- 
resented as functions of </>: V(<j>) ^ V{<f>). This simple 
method allows one to construct cosmological models de- 
scribing a smooth transition from ordinary dark energy 
to the phantom one (0, Q). 

Therefore the Friedmann equations admits the lineariz- 
ing substitution and can be studied via different powerful 



For example: 

p + 'ip 



U l = - 1 - 



U 2 = p-3p, U 3 =-(p-p), 



or 
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3p 

-2{ p + Y 



U-! = 



hp + 3p 



and so on. 

Remark 1. If the universe is filled with scalar field <j) 
whose Lagrangian is |3]) then the expression ((TJ will be 

U n = n{n - S)K + n 2 V. 

In particular case n — 3 Uz = 9V(<j>) (see ©). This 
particular case has been extensively studied in[2| , Q ■ 
Remark 2. For small values of n <C 1 one gets U n ~ 
— 3n(p + p)/2; for example, if n = 0.01 then 



-(0.0149p+0.015p) - -0.015(p + j>). 



Therefore one can use U n < to check whether the weak 
energy condition is violated [l7|. If, on the contrary, 
n>l then U n ~ n 2 p. 

In the case of general position, the solution of the equa- 
tion ([6]) has the form 

= Ci1p n + C 2 1p n , (8) 

where ip n is linearly independent counterpart of tp n : 



df 
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(9) 



Equation ([8|) is enough to establish the following theo- 
rem: 

Linearization Theorem. Let a = a(t) be the solution 
of © with k = and with p and p, given by |T|) . Then 
the two-parameter function a n — a n (t; ci, Ca): 



a n = a ci + c 2 



dt 
7577 



l/n 



(10) 
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will be solution of ([2]) with new energy density p n and 
pressure p n satisfying: 
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Another way to formulate this theorem is to say that the 
expression 



Un 



n ((n — l)d 2 + ad) 



is invariant with respect to transformation a — > a n with 
a n defined by (jTUJ) . We'll use the term "dressing" for 
the process of transformation of a triple {a,p,p}, with 
the resulting triple {a„, p n ,p n } being referred to as the 
dressed one. 

Remark 3. This theorem is valid for the case k = 0. If 
k = ±1 then this theorem will hold if and only if n = 0, 1. 
Remark 4- If n = 3/2 then the equation results in 
p n =P (but p n ^ p). 

Let a(0) = 0, i.e. suppose that at £ = there exist an 
initial singularity. Lets assume that a(£) ~ £ A for £ — ► 0. 
One might easily verify that if 2n\ < 1 then a„(0) = 0. 
Now let us choose |cl/c 2 | ^> |£max| where £ = £(£) is the 
quantity from (0) ; £ is a bounded function at the interval 
< £ < T and £ max is the maximal value of £(£) at this 
interval. It can been seen that at a given time interval 
a n (t) behaves similar to a(£) with any given precision 
rate. 

Therefore, the conjecture of our immortal astronomer 
was wrong: Her observations might result in a(£), but 
she can't be sure that the "real" scale factor is ait). It 
can be a„(£) as well. 



III. SMOOTH DYNAMICAL 
(DE)-PHANTOMIZATION 

Let us now impose the equation of state: p = wp = 
(7 — l)p, where 7 is adiabatic index. In this case the 
solution of equations @ with k — has the form 



a = a l t 2 '^\ p = 



9 7 2 £ 2 ' 



where dj =const. Using "dressing procedure" a — > a n 
one gets a new solution a n — a n (t; c\, c 2 ) such that 



t (£; ci, 0) = c\ /n a{t), o„(£; 0, c 2 ) = d^ 2 ^, 



where 



or 
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Phantom energy always takes place for negative 7 (or 
7). If, on the other hand, 7 > then a„(£;0, c 2 ) will 
describe the phantom cosmological model whenever n lies 
in the range n < or n > 37/2. For example, if initial 
a(£) describes the dust universe with w = (7 = 1) 
then a„(£;0, c 2 ) will describe the phantom cosmological 
model for n < and n > 3/2; if the initial scale factor- 
describes radiation universe with w = 1/3 then one gets 
the phantom model both for negative n and n > 2. 

If ci 7^ and c 2 7^ then (fT0)l will describe the 
smooth dynamical (de)-phantomization (see 0] for the 
case n = 3). Phantom energy results in big rip singu- 
larity at some time £ = £* . Controlling the zeroes of the 
starting solutions a(£), we can construct exact cosmolo- 
gies with the presence of big rip. In fact, let a(£*) = 
and a(t) — (£ - £*) A or a(t) — (£* - £) A for £ — £*. Us- 
ing (|10[) one can conclude that, for nX > 1, a„(£;ci,c 2 ) 
will describe the big rip as the singularuty of a(t) is ap- 
proached. 

If £(£) — > as £ — > ±00 then one has a solution allowing 
for finite "phantom regions". In other words, there exist 
two instants £1 and £ 2 , £1 < £* < £ 2 such that (p n + 
Pn)(h) = (p n +Pn)(h) = and p n +p n < if £ x < £ < £ 2 . 

All these regions can be cut out, with the edges being 
newly sewn up in such a way that neither the scale factor 
(and its first two derivatives) nor density or pressure will 
experience a jump anymore. To show this let us assume 
that the universe is born at some time instant £ = < t\ 
and suppose there exist a time interval, during which all 
the energy conditions are satisfied. One can demonstrate 
that as time goes to t%, the strong energy condition is 
the one bounded to be violated first, the weak one being 
extant. As a result, the universe would get inflated up 
to the instant t\. Furthermore, when £ = t% we arrive 
precisely to the De Sitter's phase, but this time instan- 
taneously. What follows next is the "phantom region", 
which can be cut out by matching the solutions cor- 
responding to £ = t\ and £ = £ 2 . Note that gluing the 
solutions in such a fasion we shall make sure to use the 
constraints much stricter than those usually imposed for 
the similar problems of non-relativistic quantum theory 
based on the Schrodinger equation. In fact, what the dis- 
cussed problem calls for is a C 2 gluing; namely the func- 
tions has to be matched up to their second derivatives. 
The latter requirement serves to eliminate the possibil- 
ity of discontinuities of pressure and density. A general 
solution of equation ^ is characterized by two integra- 
tion constants, and it is not obvious that three matching 
conditions will be satisfied in general. However, it is easy 
to verify that the third condition, namely the equality of 
the second derivatives, is automatically fulfilled, pro- 
vided that the first two, namely the C 1 matching has 
been successfully applied. This apparent fluke is due to 
the fact that the matching is done at an inflection point 
of log tp- 

Indeed, let us return to |5]). It is possible to write 



w 



3(tu + l) - 2n 



ipn + c^'ipn where \I> T 



) at £ < ti and 
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= v I'ri + ' > at t > t%. Using the conditions 



A. Simple generalization of k — + 1 dust model 



one can express the constants qj via c\ 2 , a(ti_ 2 ) and 
a(ti, 2 )- Since p n — n 2 d 2 /a 2 , it follows that p n (ti) = 
p n (t 2 ). But p n (t 1:2 ) = -Pn{h,2), therefore p n (*i) = 
Pnfa)- Using (JSJ) one concludes that a n (t\) = a n (t 2 ), 
a n (h) = a n (t 2 ) and a n (ti) = a„(t 2 ), and hence = 

Therefore, the finite phantom region can be cut out. 
However, an immediate question arises here: is it really 
required to do this? Of course, the smooth phantomiza- 
tion looks extremely suspicious. In the universe filled 
with scalar field </>„ (<j) n is the dressed field: <p —> <Pn) it 
leads to the highly unusual situation: <\> 2 n > at t < t\, 
t > t 2 while 2 < at < ti < t < t 2 . Is it possible for 
kinetic term to change its sign? While absolutely impos- 
sible in Minkowski space, this may still be the case in 
the General Relativity [1] (see also @,@]). Our point of 
view is this: one shall treat the smooth phantomization 
as one real physical phenomenon. If, however, the finite 
phantom region is an undesirable one then it can simply 
be cut out by the method described above. 



SPACE-TIME WITHOUT EVENTS 
HORIZON 



IV. 



The equation (|T0|) allows one to construct the space- 
time without events horizon (see also Q). In order to 
show this let tf be an instant, such that a n (tf, c\, c 2 ) = 
but a(tf) = a,f 7^ 0. In other words, suppose that 
ci + c 2 £/ = 0, where £/ = £(£/). Then, for t — > tf (and 
t < tf) one gets 



l/r, 



a„ - a f (ci + c 2 £f - c 2 (t f - t)£ f ) = — (t f - t) 

^ (12) 
where c 2 = — k 2 < 0. Integrating this equation for future 
directed radial null geodesies, ds 2 — dt 2 — a^dx 2 one will 
get 



Ax. 



aj_ 

K 



dt' 



na f 



(tf-t) 1 /" K(l-n)*^ t/ 



lim (tf-t) 



(n-l)/n 



It is easy to see that for < n < 1 we will have Ax = +00 
which shows that radial null geodesies circumnavigate 
the universe infinite number of times as the future c- 
boundary at t = tf is approached. By homogeneity and 
isotropy, we can conclude that all future endless timclike 
curves define the same c-boundary point. 

In the case n = 1 one gets A% ~ — log (tf — t) — > 
+00 as the final singularity is approached. This case is 
extremely interesting because (see Remark 3) when n = 1 
one can use the dressing procedure to construct the exact 
solutions for the universes with k = ±1. Let's consider a 
few examples. 



In the simplest dust case with p — one can solve the 
system to get 



a = a m sm rj, 2r] — smzrj — — , 

U"m. 



(13) 
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Using (Tl0|) for the case n = 1 one will obtain the gen- 
eral solution 



a{t)gen = cia(t) + c 2 a(t), 



(14) 



where c\ i2 are the arbitrary constants. It is possible to 
rewrite (|14|) in the form 

a(t) = a(t) gen = A sin 77 sin (5 - rj) , (15) 

with two arbitrary constants A and i5. This solution de- 
scribes the universe being born from the initial singular- 
ity (f?j = 0, ti = 0) and collapsing thereafter into the final 
singularity at 77/ = S or 

t f = ^-(28 -sin 25) 



1 9 6 
sin - ; 

2a m 2 



where we have introduced a new parameter a such that 

2?7 — sin2?7 = 4crf, a 
hence, 

a = . 2 Zicx sin 7? sin ((5 -77), (16) 
sin (0/2) 

and the maximum value of a = a m will occur at 77 = 5/2. 

It is easy to see that upon the choice 5 = it one gets 
the well known "dust solution" (fT5|) . In case of general 
position one shall choose < 5 < it. It can be seen that 
for7j<c/j/, 5 = 7r — e and e 1 , IT6|) will behave similar 
to (fl3|) . But if t ~ tf then one gets something really 
different: a universe without events horizon. To show 
this lets consider 

ds 2 = dt 2 — — j sin 2 ?7sin 2 (<5 — 77) [dx 2 + sin 2 x<if2 2 ] . 

Upon integration of the equation ds 2 — (describing the 
future directed radial null geodesies) one gets: 



Ax = 2 



sin 77 



sin(£-77) 
We note that if 5 — it then 



drj = +00. 



(17) 



Ax = 2(tt — rj) < 00. 

The result (fTT)) shows that radial null geodesies circum- 
navigate the universe an infinite number of times as 



5 



t — > t/. This fact and the homogeneity+isotropy results 
in conclusion that (i) this universe has no event horizons 
and (ii) all future endless timelike curves define the same 
c-boundary point. At last, 



Ax = 2 



sm rj 



sin(<5 — 77) 



drj < +00. 



One can show that all energy conditions are valid. For 
this let us point out that sum p + 2>p in general model 
is equal to the sum p + 3p in starting model (see (fTT|) 
for the case n = 1). This fact results in validity of strong 
energy condition for our model. Finally, from Friedmann 
equations one can see that density of energy is always 
positive at k = 1. By this property and the validity of 
the strong energy condition, the weak energy one will be 
satisfied automatically. 



When 8=1 the resulting solution might be denoted 
as " quasisingular" because 



lim a 

t — >oo 



lim sinh 1/2 6>lncoth- ~ exp(-0/2). 



From this relation one can see that scale factor tends to 
singularity but never achieves it. 

Both singular and quasisingular models contains no 
events horizons. To show this for the quasisingular case 
let us integrate the equation for future directed radial 
null geodesies (ds 2 = 0) just like it has been done in the 
previous subsection: 



Ax 



d6 



sinh 1/2 (91ncoth£' 



B. Generalization of a Lambda-radiation model in 
flat space 

Let's consider the flat universe which has a positive 
vacuum energy A. Let's also assume that the universe 
is filled with the radiation. Solving system @ we will 
obtain the initial solution for the scale factor 



a = do s 



sinh 1/2 6, = 2VXt, 



p = A(l + sinh" 



(18) 



where ao is a positive constant. If t > t v = 
^/^arcsinh(l), the strong energy condition will neces- 
sarily be violated. Using ([10]) for the case n = 1 one can 
see that general solution can be written in the following 
form 



1/2 coth 6 | 
a = ao sinh ' 6 In — - . 



(19) 



where 8 is a positive constant. The parameter e, intro- 
duced here, plays an important role in our reasonings. 
If e = 0, then (fTO)) will be equivalent to the initial so- 
lution. In the remaining cases one can without any loss 
of generality assume e = 1. There will be three types of 
solutions. If 8 < 1 the universe will be open. This type 
of solutions has the following asymptotic behavior 



a — > — 0.5ao In 8 exp(v At), t —* 00. 

It is easy to see that this universe is plagued by the events 
horizon. 

The case 8 > 1 is a more interesting one. If 8 « 1, 
then solution will describe the universe, starting from an 
initial singularity (0 = 0, t = 0) and ending up in the final 

singularity at tf = y^arccoth<5. One shall note that, if 

8 > 1 + v2 the strong energy condition will always be 
satisfied (in fact, universe will end up in singularity long 
before the time t„). 



This integral diverges because subintegral expression has 
an exponential asymptot at large 9. Therefore radial null 
geodesies circumnavigate the universe an infinite num- 
ber of times as t — > 00. This fact and the homogene- 
ity+isotropy result in the conclusion similar to the one 
from the subsection B, namely that such universe pos- 
sess no events horizon. Absence of events horizon for 
singularity model can be proved by analogy. 

In conclusion let us analyze the equation of state for 
the generalization of a lambda-radiation model. One can 
show that the value w = p/p is equal to 



w 



1 2 (1 -sinh 2 6»)ln 



2 coth ' 



3 3 (cosh 8 In 



coth - 



2) 2 



(20) 



From this relation it follows immediately that for both 
open and quasisingular models w — * —1 at large t (large 
6). For the closed model w = —1/3 whenever we ap- 
proach the final singularity. In the initial singularity 
W = 1/3 for all models. Close examination of equation 
([2H)l shows that w is always greater than -1 for all cases, 
i.e. weak energy condition will always be satisfied. 



V. POSSIBLE GENERALIZATIONS OF THE 
METHOD 

In the previous sections we have shown that no finite- 
time observations made by our astronomer will yield an 
exact picture of the future dynamics of her universe. One 
can claim, however, that the cited " nonuniqueness" is 
due to the special properties of the Friedmann equations 
and will probably disappear for the cosmological models, 
distinct from this one. In current section we are going 
to show that, in fact, the above-mentioned ideology can 
be applied absolutely similarly in at least two particular 
cases: for anisotropic universes and for a certain brane 
world models. 
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A. Anisotropic universes 
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c 




-1 


3 




c- 




[i- 




2 



In this subsection we are going to show that the de- 
scribed method can be generalized to some cases of 
anisotropic universes. Consider a simplest possible case 
of anisotropic model: 



ds 1 



dt 2 - a 2 (t)dx\ - b 2 (t)dx\ - c 2 {t)dx\, 



with the perfect fluid energy-momentum tensor. The cor- 
responding system of equation has the form 




and 



a 
a 



(P-P), 



\{P-P) 



(p + 3p). 



(21) 



Using (|2T|) one can rewrite (|22|) as 

dbc + deb + abc — Spabc. 



(22) 



(23) 



There are known quite a few exact solutions of the system 
(|2"Tj) . (f2"2"|) . such as the vacuum solution of Kasner 
or dust solution of Safiicking, Heckmann Our aim, 
however, is to develop the dressing procedure for the (|2Tj) . 

In order to do this lets first define a new function 
ip = abc. Using (f2"Tj) one can check this function to be a 
solution of the Schrodinger equation: 



4> = Vi/i, 



(24) 



with potential V — 9(p — p)/2. In fact, this is just a 
replica of the case n = 3 of the Friedmann equation. 
The general solution of the equation (|24|) is 



(25) 



where ip is linearly independent counterpart of -0. Since 
^ is the solution of the with initial potential V (with 
the new density and pressure p and p, but same differ- 
ence: p — p = p — p), we can define the new solutions of 
the system ([2~Tj) A, B, C with similar right hand side: 



(26) 



A 


A 


' B 
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~A ^ 


-A 




■s) 


~~ 2 


B 


B 
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B 


-B 


(4- 




2 



ip-p), 



(27) 



and 



(P-P), 



ABC + ACB + ABC = 3/5*. 



(28) 



(29) 



Further on, one can represent equations (|26|) , (|27|) as the 
particular system 

2AAV + 2AA^ -2^ A 2 +3(p- p)A 2 * = 0, , , 
2BB* + 2 J BB*-2*B 2 + 3(p-/7)B 2 * = 0. ^ ' 

One can show that if A, B and * are solutions of respec- 
tively (M and mj then the function 



C = 



(31) 



would be a solution of (|28[) . Thus, assuming p, p and 
* to have certain values, and solving (|30[) for A and -B, 
one can automatically calculate C (via (|3"Tj) ). The new 
density p will then be defined from the eq. (f29|) and the 
new pressure from 

P = P- P + P- 

Equations (|30[) can be solved by simple substitution A = 
exp(a), B = exp(/3) which results in system 

* 3 

* • 3 



therefore 



and 



A = F(t) exp ^c, J 



1 dt 1 



B = F(t)ez P (c II f^ 



where 

F(t) = exp 



dt' 



dt" (p(t")-p(t")) 



2 J *(t' 
c r and c„ are arbitrary constants. 

B. A Few notes about the brane universe 



Let's start our considerations from the following sys- 
tem: 

2 



2XJ 



-2-=p + 3p+^(2p + 3p). 
a A 



(32) 
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where p is the density, p 
on the brane. 



pressure, and A is the tension a(Q) = then 



Before we actually start applying our method, let us 
make one simple but interesting observation. If we'll try 
to solve the system (|3"2"j) for the p as a function of a (for 
finite values of A ^ 0), we will unexpectedly end up with 
not just one, but two distinct solutions: 



pi 2 = -A± WA 2 + 2A 



(33) 



For the negative values of A this (interestingly enough) 
gives an upper bound on the absolute value of the Hubble 

constant: H 2 = (f J < —A/2 and presents two positive 
values of p. In the case of positive A (which in theory 
might correspond to our universe) one of the values of p 
will be negative (it might seem strange on a first glance, 
but the study of the models with the negative p has expe- 
rienced quite a surge in a past few years. See, for example 
[l(|). In other words, even the exact knowledge of the 
very functional representation of scale factor would not 
allow the astronomer (living in such universe) to uniquely 
determine the density function! 

With this being said let us finally state the following 
version of our linearization theorem: 

Linearization Theorem for Branes. Let a = a(t) 
be the solution of (|32| with k = and with p and p being 
the density and pressure functions correspondingly. Then 

1. Function ip n = a n will be solution of the Schrodinger 
equation 



a(t) = a(to) 



H 



i(3A 7 i + 2V2A) 



1/37 



\t (3\~ft Q + 2 V / 2A) / 



P = 



a 3 7 i(3VA 7 i + 2V2) ' 
4A 

37i(3\/A7t + 2\/2) ' 



(37) 



For the dust case w = (7 = 1) and n = 3 one gets the 
following two-parameter function 03 = 03 (t; ci, ci): 



Cl 



2«(ffil + 3i 

_ (38) " 

where ((t) = 3\t + 2\/2A. One can show that in this 
case the solution might above all describe the space-times 
without events horizon (for a special choice of values of 
ci, c 2 ). 

Remark 5. If 7 < then one get the phantom model with 
two big rips jlQ ]: 



1/3 



a(t) = 



/'" 



1 -4.5A7 2 (t-T) 5 



l/3|7l 



where p? and T are constants. Using our method one can 
obtain two-parameter (or 2iV-parameter) function a n — 
a n (t; ci, C2) such that a n (t; c\ = 1, C2 = 0) = ait). 



with the potential 

W n = I (2np-3(p + p) + ^(np-3(p + p))) 



(34) 



1 



2n{K + V) - 6K + -{K + V) (n(K + V) - 6K) 
A 



(35) 

(Note, that despite the discussed non-uniqueness of the 
density function, potential W n is unique!) 

2. The two-parameter function a n = a n (t; c\, c 2 ), satis- 
fying (flC)]) will be solution of (|32|) with new energy density 
p n and pressure p n satisfying: 



W n (p,p) = W n (pn,Pn)- 



(36) 



As before, the proof is nothing more but a technical 
exercise that will be left out for our readers. 

Example. System (|3"2"j) can easily be solved when 
p/p = w =const. Let A > 0, w > — 1 (7 > 0) and 



VI. CONCLUSION 

In this paper we have discussed a simple (and easily 
automatizable) method of construction of exact solutions 
of Friedmann equations. Despite simplicity, the method 
allows for acquirement of solutions characterized by ex- 
tremely interesting properties. What is more, it appears 
(for k=0) that the very abundance of the set of solutions 
that are to be obtained this way leads us to a stunning 
conclusion: no matter how accurate our astronomical ob- 
servations are, there exist not just one, but a whole set 
of solutions that will satisfy the observational data while 
leading to essentially different dynamics in future. This 
sudden twist leads us to seemingly unavoidable conclu- 
sion about the principle indefiniteness of the future, hid- 
den in the Fridmann equations. For a first glance such 
conclusion looks really disappointing, rendering useless 
all our efforts to build a suitable cosmological model de- 
scribing our universe. 

However, everything above-said doesn't mean the im- 
possibility to determine the actual dynamics of the uni- 
verse in principle. Even though the usual observational 
methods might not give us the final answer, we can use 
the statistical methods to get it - that is, anthropic prin- 
ciple. 
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For example, in the article we have studied the partic- 
ular models possessing the final singularity and no event 
horizons. What is the practical significance of such mod- 



els? As follows from the recent observations (see 11 1, 
12 1 ) our universe suffers the accelerated expansion 13| , 
14j . As for now, the most probable cause of such expan- 
sion lies in nonzero cosmological constant. If this is really 
the case, then the future dynamics of observable universe 
is confined within the particles horizon and, as such, leads 
to problems with formulation of a fundamental physical 
theory (like the string theory or hypothetical M-theory) 
in a finite volume [l5| . If, on the other hand, the lifetime 
of the universe (with an observable expansion rate) will 
exceed the limit of 10 60 years, the dominating observers 
will, as follows from [l6| , be the ones of a quantum fluc- 
tuations origin, which, of course, could hardly be called 
compatible with our observations. Therefore, as follows 
from the anthropic reasons of [l6j], the most probable 
scenarios would be those describing the contemporary 
expansion, being traced by the consequent contraction 
phase and the "horizonless" collapse [19| . In other words, 
exactly that type of scenarios, that has been naturally 
constructed with the aid of "dressed" initial solutions. 



Finally, let us emphasize that the purpose of this paper 
was a mere study of Friedmann equation's solutions that 
has been constructed in this article by the linearization 
procedure. The practical aspects of application of the 
described technique to the observed universe lied out of 
scoop of this work. However, we are hoping to get back 
to it (as well, as to applications of the method to more 
interesting anisotropic and brane models) in one of our 
future articles. 
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